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ABSTRACT 



We report the discovery of a bow shock around the high-mass X-ray binary (HMXB) 4U 1907+09 using the Spitzer Space Telescope 
24 pm data (after VelaX-1 the second example of bow shocks associated with HMXBs). The detection of the bow shock implies that 
4U 1907+09 is moving through the space with a high (supersonic) peculiar velocity. To confirm the runaway nature of 4U 1907+09, 
we measured its proper motion, which for an adopted distance to the system of 4 kpc corresponds to a peculiar transverse velocity of 
=! 160 +115 kms -1 , meaning that 4U 1907+09 is indeed a runaway system and supporting the general belief that most of HMXBs 
possess high space velocities. The direction of motion of 4U 1907+09 inferred from the proper motion measurement is consistent with 
the orientation of the symmetry axis of the bow shock, and shows that the HMXB is running away from the Galactic plane. We also 
present the Spitzer images of the bow shock around VelaX-1 (a system similar to 4U 1907+09) and compare it with the bow shock 
generated by 4U 1907+09. 
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1. Introduction 
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4U 1907+09 is a high-mass X-ray binary (HMXB), consisting 
of an X-ray pulsar (neutron star) accreting from a late O-type 
supergiant companion star (Giacconi et al. 119711 Schwartz et 
al. [TM)l Marshall & Ricketts [19801 Maki shima et al. [19841 
van Kerkwijk, van Oijen & van den Heuvel [T989t Cox, Kaper 
& Mokiem 125531 Nespoli, Fabregat & Mennickent l2558l ). Like 
other HMXBs, 4U 1907+09 is the descendant of a massive bi- 
nary, whose primary (initially more massive) star experienced 
supernova explosion (van den Heuvel & Heise[1972). The fact 
that the binary has remained bound after the supernova explo- 
sion implies either that the system lost less than a half of its 
initial (pre-supernova) mass (in the case of a symmetric super- 
nova explosion; Boersma ll96H or that the kick attained by the 
neutron star at birth had a favorauble orientation and magnitude 
(if the explosion was asymmetric; Hills 119831 Tauris & Takens 
11998b . In both cases, the HMXB attains a recoil velocity, which 
could be large enough (i.e. > 30 kms -1 ) to classify the sys- 
tem as a runaway (Blaauw 1 19611 1. The distribution of HMXBs 
in the Galaxy and their kinematic properties are consistent with 
the possibility that most of them are runaways (van Oijen [T9891 
Chevalier & Ilovaisky |1998l >. A similar conclusion was derived 
by Coe d2"5531 > for HMXBs in the Small Magellanic Cloud. 

The high peculiar velocities of HMXBs could be revealed 
directly, via measurement of their radial velocities (van Oijen 
1989) and/or proper motions (Chevalier & Ilovaisky 1998), or 
indirectly, through the detection of the secondary attributes of 
runaway stars - the bow shocks (e.g. van Buren & McCray 
1988; van Buren, Noriega-Crespo & Dgani 1995; Gvaramadze 



& Bomans 2008; Gvaramadze, Kroupa & Pflamm-Altenburg 
2010a; Gvaramadze, Pflamm-Altenburg & Kroupa 2011). The 
first-ever detection of a bow shock associated with a HMXB was 
reported by Kaper et al. (1997), who discovered a very sym- 
metric Ha arc around the HMXB VelaX-1 (a system similar 
to 4U 1907+09). This discovery unambiguously showed that at 
least in the case of VelaX-1 the system got its high space ve- 
locity owing to the binary-supernova explosion [an alternative 
mechanism for the origin of runaway stars (including binaries) 
is based on dynamical three- and four-body encounters in dense 
star clusters; Poveda, Ruiz & Allen [19671 Gies & Bolton [T986l >l. 

Huthoff & Kaper ( 120021 ) used the high-resolution IRAS maps 
to search for bow shocks around 1 1 high- velocity HMXBs and 
detected only an infrared counterpart to the already known 
bow shock generated by Vela X- 1 . The paucity of bow shock- 
producing HMXBs is consistent with the observational fact that 
only a minority of runaway OB stars are associated with (de- 
tectable) bow shocks (van Buren et al. 119951 . The most reliable 
explanation of this fact is that the majority of runaway stars are 
moving through a low density, hot medium, so that the emis- 
sion measure of their bow shocks is below the detection limit 
or the bow shocks cannot be formed at all because the sound 
speed in the local interstellar medium (e.g. the hot gas within 
large-scale bubbles around OB associations) is higher than the 
stellar space velocity (Kaper, Comeron & Barziv [T999l Huthoff 
& Kaper[2002l. In case of bow shocks with small emission mea- 
sure, one might expect to detect them using modern infrared tele- 
scopes. 

Motivated by the above arguments, we have undertaken a 
search for bow shocks around the HMXBs from the sample of 



1 



V.V. Gvaramadze et al.: 4U 1907+09: a HMXB running away from the Galactic plane 



Huthoff & Kaper ( 120021 see their Table 1) using the archival 
Spitzer Space Telescope data. We have discovered a bow shock 
around the HMXB 4U 1907+09 and detected the already known 
bow shock associated with VelaX-1 (Sect. [3]). The orientation of 
the bow shock generated by 4U 1907+09 is consistent with ejec- 
tion of this HMXB from the Galactic plane. To confirm the run- 
away nature of the system, we have measured its proper motion 
using archival data over a total time baseline of 50 yr (Sect.|4|. 
Sect. [5] discusses the results obtained in the previous sections. 
We summarize and conclude in Sect. [6] 

2. 4U 1907+09: general data 

The HMXB 4U 1907+09 was disc overed with UHURU as a dis- 
crete X-ray source (Giacconi et al. 11971) . The X-ray source was 
subsequently identified with a highly reddened early-type star 
with broad, strong Ha emission (Schwartz et al. 1980). Schwartz 
et al. d 19801 ) hypothesized that 4U 1907+09 is a binary system 
composed of an OB supergiant star and a compact object, and 
that accretion of the stellar wind onto the compact object is re- 
sponsible for the X-ray emission. This hypothesis was proved 
by Marshall & Ricketts ( 1980), who revealed a binary period of 
- 8.38 day in the X-ray light curve of 4U 1907+09, confirm- 
ing that the system is a new member of a class of HMXBs (see, 
e.g., Kaper & van der Meer l2007l for a review on these objects). 
Makishima et al. ( 119841) discovered a 437.5 s X-ray pulsations 
from 4U 1907+09, establishing this system as a binary X-ray 
pulsar. 

Follow-up optical spectroscopy of 4U 1907+09 by Cox et 
al. (|2005 ) showed that the companion of the X-ray pulsar is a 
late O-type supergiant star (08 Ia-09 la; cf. van Kerkwijk et al. 
119891 lye 1 1986b . This classification was refined by Nespoli et 
al. (120081 ). who got the infrared spectrum of 4U 1907+09 and 
derived the spectral type of 09.5 lab. 

The distance to 4U 1907+09 still remains uncertain (cf. van 
Kerkwijk et al. [19891 Cox et al. 120051 Nespoli et al. 12008b . 

Using the optical photometry from Schwarts et al. (11980) . 
V = 16.37 + 0.02 mag, B - V = 3.17 + 0.06 mag, and adopt- 
ing the intrinsic (B - V)o colour of an 09.5 1 star of -0.26 mag 
(Martins & Plez 120061 ). one finds the interstellar reddening to- 
wards 4U 1907+09 E(B - V) = 3.43 + 0.06 mag. Then taking 
Ay = 3AE(B - V), one obtains A v 1 0.6 + 0.2, which for 
M v = -6.34 + 0.45 mag (Martins & Plez l2006l leads to a dis- 
tance d =: 2.6^q j kpc. 

Similarly, using J and K s magnitudes of 4U 1907+09 
(10.00 + 0.02 a nd 8 7 7 + 0.02, respectively) from 2MASS 
(Skrutskie et al. I20061 . (J - K s ) = -0 .21 m ag and M Ks = 
-5.52 ± 0.20 mag from Martins & Plez (12006b, an d adopting 
the extinction law from Rieke & Rebofsky (1985), one finds 
A Ki = 0.66[(7 - K s ) -(J- KM = 0.95 + 0.0^1 suggesting 
a distance of =* 4.7^H kpc. The large discrepancy between the 
distance estimates based on the visual and near infrared photom- 
etry could be understood if the reddening towards 4U 1907+09 
is anomalous, i.e. the total-to-selective absorption ratio Ry < 3.1 
(cf. Cox et al. 120051 . 

Using a similar procedure and adopting the intrinsic colours 
and absolute magnitudes from Wegner (119941 12006) . Nespoli 



1 We also tried to estimate A Ks using the Spitzer photom- 
etry of 4U 1907+09 and the Stellar Performance Estimation 
Tool ( http://ssc.spitzer.caltech.edu/warmmission/propkit/pet/starpet/ in- 
dex.html), but found that the fluxes predicted by this tool (based on 
Kurucz-Lejeune atmospheric models) cannot be reconcilled with the 
observed ones. 



et al. d20TJ8b derived a distance to 4U 1907+09 of 2.8^° kpc, 
where the quoted errors stem mainly from the uncertainty in the 
absolute magnitudes. 

A distance d = 2.1 - 2.6 kpc was derived by van Kerkwijk 
et al. d 19891) and Cox et al. (12005) from interstellar absorption of 
sodium in the spectrum of 4U 1907+09, while the detection of 
a high velocity (~ 60 kms~') component in the Ki absorption 
line in the line-of-sight towards 4U 1907+09 suggests a distance 
of > 4 kpc (Cox et al. 12005) . The latter estimate is based on 
the model of Galactic rotation by Brand & Blitz ( 11993) . which 
assumes that, to first order, the Galactic velocity field is axisym- 
metric. Note, however, that 4U 1907+09 is located in the direc- 
tion of the Sagittarius arm tangent region and that our line-of- 
sight towards this HMXB passes through the spiral arm over a 
distance of several kpc. One cannot, therefore, exclude that ran- 
dom (non-circular) motions (up to several tens of km s ) caused 
by outflows from star-forming regions could contribute to the 
origin of the high-velocity components in the K i absorption line, 
thereby affecting the kinematic distance estimate. 

Finally, we note that the interstellar extinction cannot be con- 
strained by using the empirical relationships between Ay and the 
interstellar column density Nu (derived from modelling the X- 
ray spectrum of the neutron star), because the total column den- 
sity towards 4U 1907+09 contains a variable contribution from 
the circumstellar material of a magnitude comparable to that of 
the interstellar one (Cook & Page 1987 in 't Zand, Strohmayer 
& Bavkal [1997l Rivers et al. 1201731 . 

As a compromise between the different distance estimates, 
we adopt a distance to 4U 1907+09 of 4 kpc in the following. 

3. 4U 1907+09: bow shock 

To search for bow shocks around HMXBs from the sample of 
Huthoff & Kaper (2Q02I, we utilized the imaging data from the 
Spitzer Space Telescope archive, obtained with the Multiband 
Imaging Photometer for Spitzer (MIPS; Rieke et al. 12004b and 
the Infrared Array Camera (IRAC; Fazio et al. 12004) . The 
publicly available data cover the fields containing 5 of the 1 1 
HMXSs (namely, GX301-2, 4U 1907+09, CygX-1, VelaX-1 
and V615 Cas). Visual inspection of the MIPS and IRAC images 
results in a discovery of a bow shock generated by 4U 1907+09. 

FigurepJ shows the MIPS 24 pm, IRAC 8, 5.8, 4.5 and 3.6 
pm, and DSS-II (red band) imagefl of the field around the 
HMXB 4U 1907+09. Like many o ther bo w shock s discovered 
with Spitzer (e.g. Gvaramadze et al. 120 lOal b. 1201 II Gvaramadze 
& Gualandris l201 lb . the bow shock associated with 4U 1907+09 
is only visible at 24pm (although the IRAC 5.8 and 8 pm im- 
ages show the gleam of emission possibly associated with the 
bow shock). The bow shock has a clear arcuate shape with the 
apex at 0'J from 4U 1907+09, which at the distance to the 
system of 4 kpc corresponds to the linear separation R \, s 0.8 
pc (see Sect. 15.3) . 

We also detected the already known bow shock associated 
with V elaX-1 [see Kaper et al. ( 119971 ) and Huthoff & Kaper 
(120021) for the Ha and IRAS 60 pm images of the bow shock, 
respectively]. Unlike 4U 1907+09, the bow shock produced by 
VelaX-1 is visible not only at 24 pm, but also at all four IRAC 
bands, i.e. at 3.6, 4.5, 5.8 and 8 pm. In Fig.[2jwe present for the 



2 The MIPS and IRAC images were obtained within the framework 
of the 24 and 70 Micron Survey of the Inner Galactic Disk with MIPS 
(MIPSGAL; Carey et al. 12009) and the Galactic Legac y Infrar ed Mid- 
Plane Survey Extraordinaire (GLIMPSE; Benjamin et al. 12003) . respec- 
tively. 
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Fig. 1. From left to right, and from top to bottom: Spitzer MIPS 24 pm, IRAC 8, 5.8, 4.5 and 3.6 pm, and DSS-II (red band) images 
of the field around the HMXB 4U 1907+09. The arrow shows the direction of peculiar velocity of 4U 1907+09, as suggested by the 
proper motion measurement (see text for details). A bright infrared source to the north is the maser source IRAC 19071+0946. In 
all panels the position of 4U 1907+09 is indicated by a circle. 



E N 


•> *• - * . 




* 






v 


... V ' 


2 arcmin 





Fig. 2. Spitzer MIPS 24 pm 
(left) and IRAC 8yum images 
of the bow shock around the 
HMXB VelaX-1. The arrow 
shows the direction of pecu- 
liar velocity of VelaX-1 (see 
Sect.|53]for details). 



first time the MIPS 24 pm and IRAC 8 pm images of the bow 
shock (Program Id.: 30088, PI: A.Noriega-Crespoj^ showing its 
fine structure (most prominent at 8 pm) similar to that detected 
in Ha by Kaper et al. ( 119971 ). Using the IRAC 8 pm image of 



3 Note that the bow shock produced by Vela X- 1 was independently 
observed with Spitzer by R.Iping (Program Id.: 30174). 



the bow shock, we estimated R \, s ^ 0.5 di$ pc, where d\,g is the 
distance to VelaX-1 in units of 1.9 kpc (Sadakane et al. |19851 l. 

The symmetry axis of a bow shock generated by a super- 
sonically moving star reflects the direction of motion of the star 
with respect to the ambient medium, which also could be in 
motion relative to the local standard of rest at the location of 
the star. For example, high-velocity outflows driven from young 
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Galactic longitude 

Fig. 3. 3° x 2° MSX 8.3 pm image of the Galactic plane (shown 
by a dashed line) centred at / = 43°.5, b = 0°, with the position of 
the HMXB 4U 1907+09 indicated by a circle. The arrow shows 
the direction of motion of 4U 1907+09. The orientation of this 
figure is the same as for Fig.Q] The brightest nebula on the image 
is the background star-forming region W 49A. 



star clusters by the collective effect of stellar winds can create 
bow shocks around (low-velocity) stars in the cluster's halo (e.g., 
bow shocks around LL Ori and several other stars in the Orion 
Nebula, all of which are faced towards the central Trapezium 
cluster; Bally, O'Dell & McCaughrean 2000). For runaway stars 
located far from star-forming regions and their associated gas 
outflows (caused by stellar winds and supernovae), one can as- 
sume that the peculiar velocity of the ambient medium is negli- 
gible compared with the space velocity of the star, meaning that 
the symmetry axis of the bow shock coincides well with the di- 
rection of stellar motion. In this case, the orientation of the bow 
shock can be used to back-trace the trajectory of the star to the 
parent cluster even for those stars whose proper motions cannot 
be measured with a high confidence (Gvaramadze & Bomans 
120081 Gvaramadze et al. l2010al b. E0TH . 

Figure|3] shows the Midcourse Space Experiment (MSX) 
satellite (Price et al. [200Tb image of a 3° x 2° region of the 
Galactic plane (centred at I = 43°5, b = 0°), with the position 
of 4U 1907+09 indicated by a circle (/ = 43°7436, b = 0?4753). 
The symmetry axis of the bow shock around 4U 1907+09 (see 
Fig.Q] whose orientation is the same as for Fig. [3} suggests 
that this HMXB is moving away from the Galactic plane. To 
prove unambiguously that the origin of the bow shock is due 
to the high space velocity of 4U 1907+09, we searched for ex- 
isting proper motion measurements for this system using the 
Vizier catalogue access toofl. We found three catalogues giving 
proper motions for4U 1 907+0 9, USNO-B1 (Monet et al. 120031 . 
UCAC3 (Za charias et al. EOTOt and PPMXL (Roser, Demleitner 
& Schilbach l20101 >. All three give huge proper motions, which, 
if real, would correspond to the peculiar velocity of several 
thousands of kms -1 . We discovered, however, that USNO-B1 
has two entries separated by less than 2" at the position of 
4U 1907+09 which (possibly) originate from mismatching of 
plates from different epochs and causes the huge proper mo- 
tion. PPMXL followed the cross-matching of USNO-B1, hence 
the star is one of the 10% (Roser et al. 12010) double entries in 
PPMXL. Also, it turns out that the UCAC3 proper motion is 
marked "with doubts" because it relies on "less than two good 



matches". Thus, an individual determination of the proper mo- 
tion of 4U 1907+09 became needful. 



4. 4U 1907+09: proper motion measurement 

With its intermediate optical magnitudes, 4U 1907+09 is not 
saturated on photographic Schmidt plates (see Fig.[TJi, provid- 
ing different epoch observations over five decades, and it is 
well measured in the CCD-based Carlsberg Meridian Catalog 
14 (CMC14fl However, there are several difficulties in the field 
around the target: 1) overlap of star images, changing already 
with the different optical passbands and leading to a completely 
different appearance of the field in the near infrared due to strong 
reddening, 2) absence of reference galaxies for direct absolute 
proper motion measurement, 3) different location of the target 
on the Schmidt plates resulting in different residual geometric 
distortions after any full Schmidt plate calibration. 

We collected all available Schmidt plate data, POSS-1 blue 
(O) and red (£), POSS Quick V, and POSS-2 B,R, and /, al- 
together 9 plates, and extracted small 15x15 arcmin 2 FITS im- 
ages from the DSS0 In case of the POSS-1 red (E) plate we also 
extracted the FITS fi le from the SuperCOSMOS Sky Surveyfl 
(Hambly et al. 120011 ). which we preferred because of its higher 
resolution. To overcome the above mentioned difficulties, we 
used 14 reference stars (of magnitudes similar to those of 4U 
1907+09) closely distributed around the target within about 3 
arcmin for a new local astrometric calibration of the Schmidt 
plates with respect to the CMC14 positions of the reference 
stars. The astrometric measurements and transformations from 
the Schmidt plates to the CMC14 system were carried out using 
the ESO SkyCat Toofl To avoid colour-dependent effects, we 
finally preferred to combine only the E,R and V plates with the 
CMC14(r')-band data. 

The resulting proper motion from simple linear fitting of the 
target positions at five epochs (1952.395, 1983.611, 1987.411, 
1992.586, 2001.956) is p a cosS = -4.1 + 5.5 masyr 1 ,^ = 
4.9 + 4.4 mas yr -1 . For the absolutisation of this relative proper 
motion we investigated the proper motions of the reference stars 
as given in the PPMXL catalogue (Roser et al. [2010), which is 
on the International Celestial Reference System. The individ- 
ual PPMXL proper motions of the 14 stars are all smaller than 
1 mas yr~ 1 , and their mean PPMXL proper motion is p a cos 8 = 
-3.6 + 3.0 masyr _1 ,jUa = -3.4 + 3.0 mas yr~'. The quoted errors 
are conservative estimates including the systematic PPMXL un- 
certainty of 1 to 2 mas yr~' in any small field (~1 square degree) 
of the sky due to the small number density of Hipparcos stars 
(Roser et al. l2010l l. Taking into account the mean proper motion 
of the 14 reference stars we then computed an absolute proper 
motion of 4U 1907+09: p a cosd = -7.7 + 6.3 masyr 1 ,^ = 
1.5 + 5.3 masyr -1 , or in Galactic coordinates: 



-2.2 + 5.5 mas yr , pi, — 7.5 +6.1 mas yr 



5. Discussion 

5.1. 4U 1907+09 as a runaway HMXB 

To convert the observed proper motion of 4U 1907+09 into 
the transverse peculiar velocity, we used the Galactic constants 



http://webviz.u-strasbg.fr/viz-bin/VizieR 



http ://w w w. as t. cam . ac.uk/cmt/ cmc 1 4 . html 

6 http://archive.stsci.edu/cgi-bin/dss_plate_fmder 

7 http://www-wfau.roe.ac.uk/sss/pixel.html 

8 http://archive.eso.org/cms/tools-documentation/skycat 
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R = 8.0 kpc and @ = 240 km s" 1 (Reid et al. 2009) and 
the solar peculiar mot ion (U @< Ve, W Q ) = (10.0, 11. 0,7.2) kms -1 
(McMillan & Binnev [20 101 )1. We found the components of the 
peculiar velocity: 

vi = 59 + 104 kms -1 , vj = 150+116 kms" 1 (1) 

(here we assumed that the errors are due to errors in the proper 
motion measurement). Taken at face value, these velocity com- 
ponents suggest that 4U 1907+09 is running away from the 
Galactic plane in the direction consistent with the orientation of 
the symmetry axis of the bow shock (see Figs[3]and[T] recall that 
the orientation of both figures is the same). 

By tracing back the trajectory of 4U 1907+09 and assuming 
that this HMXB originated in the Galactic plane, we inferred that 
the system was ejected from the region around I = 43°.5, b — 0°. 
We did not find, however, in this region (of radius of 1°) any 
known young (S; 10 Myr) star cluster or association located at 
a distance comparable to that of 4U 1907+09. (The only young 
star-forming region in this area is the background one W 49A. ) 
This could simply be owing to the huge interstellar extinction 
in the Galactic plane. Another possibility is that the progenitor 
massive binary was dynamically expelled from the parent cluster 
well before the primary exploded as a supernova and that the sys- 
tem changed the direction of its motion because of a kick caused 
by the supernova explosion (cf. Pfiamm-Altenburg & Kroupa 
120101 1. In this case, the parent cluster could be located at sev- 
eral degrees from the region inferred from the proper motion 
measurement and orientation of the bow shock. 

We also searched for the parent diffuse supernova remnant 
in the same area of the Galactic plane using the catalogue of 
Galactic supernova remnants by Green ((2009 0. The only su- 
pernova remnant listed in the catalogue in this area is the back- 
ground one W49B. The non-detection, however, is not surpris- 
ing because the kinematic age of 4U 1907+09 (i.e. the time 
elapsed since the supernova explosion), f^n = dsinb/vi, ~ 10 5 
yr, is comparable to the typical lifetime of diffuse supernova 
remnants (e.g. Shull, Fesen & Saken [T989l . 



Roche lobe, R 2 ~ Rl, where R^ is the radius of the Roche lobe, 
given by (Eggleton ll983l) 

0.49 apre - SN g 2 / 3 

L 0.64 2 / 3 +ln(l+<? 1 / 3 )' {> 

and q = M 2 /Mi. Eq. (0 also shows that the larger M e j = M\ - 
M co the higher v sym . Assuming that at the moment of supernova 
explosion the secondary was a 26 M main-sequence star (cf. 
Cox et al. |20051l of radius 

fl 2 = O.8(M 2 /M ) (X7 R (4) 

(Habets & Heintze 119811 ) and adopting the maximum pre- 
supernova mass of the exploding star of 17 M [see Fig. 6 
of Meynet & Maeder (120031 )1. one finds from Eqs(0-(0 that 
flpre-SN — 2.4/?l — 20 R and the maximum possible systemic 
velocity v™^ =* 230 kms -1 . The latter figure, when compared 
with the total peculiar (transverse) velocity of 4U 1907+09, Vtr = 
(v^ + v 2 ,) 1 ' 2 - 160+115 kms~' [see Eq. Q}], means that the sym- 
metric supernova explosion could be responsible for the origin 
of the high space velocity of 4U 1907+09. 

The maximum recoil velocity would be smaller if at the mo- 
ment of supernova explosion the secondary star already entered 
the supergiant stage of evolution. Adopting R^ ~ R 2 - 26 R 
(Cox et al. 120051 . one finds v™ y ™ - 130 kms -1 . Still, the "mea- 
sured" peculiar velocity of 4U 1907+09 within the l<x error bar 
is consistent with the maximum possible recoil velocity. Note, 
however, that a pre -sN is related to the post-supernova semimajor 
axis, flpost-SN, through the relationship: 

(M 2 +2M co -MA 

flpie-SN - — — Apost-SN , 

\ M co + M 2 ) 

from which it follows that a pie -SN - 24 R if one uses the stel- 
lar masses adopted above and the present semimajor axis of 
4U 1907+09 of =s 50 R . This estimate suggests that at the mo- 
ment of supernova explosion the secondary star was still on the 
main-sequence and left it later on. 



5.2. Origin of a high peculiar velocity of 4U 1907+09 
5.2.1 . Symmetric supernova explosion 

Boersma (1961) showed that a post-supernova binary remains 
bound if the mass of the supernova ejecta, M e j, comprises less 
than a half of the system's initial mass, i.e. M e j < (l/2)(Mi + 
M 2 ), where Mi and M 2 are the pre-supernova masses of the pri- 
mary and secondary stars, respectively. In response to the instant 
mass loss caused by the supernova explosion, the binary system 
recoils with a velocity, given by (Gott, Gunn & Ostriker 119701 
Iben & Tutukov [1997j l 

_/ GM 2 \ 1/2 / M 2 \' /2 M! -M co 

VSym \ Opre-SN / \Ml+M 2 ) M 2 +M co ' 

where G is the gravitational constant, a pre -sN is the pre- 
supernova semimajor axis of the binary, and M co = M\ — M e j is 
the mass of the compact object (a neutron star or black hole). In 
case of 4U 1907+09, the stellar supernova remnant is a neutron 
star and we assume that it has a canonical mass M co = 1.4 M . 

It follows from Eq. (0 that the HMXB achieves the highest 
peculiar velocity if the pre-supernova binary was as tight as pos- 
sible, i.e. if the secondary star of radius R 2 was close to filling its 

9 http://www.mrao.cam.ac.uk/surveys/snrs/ 



5.2.2. Asymmetric supernova explosion 

The recoil velocity could be higher if the supernova explosion 
was asymmetric and the stellar remnant received a kick at birth. 
In this case, the systemic velocity of the post-supernova binary 
depends also on the magnitude and orientation of the kick, and 
is given by (Stone [19 821 ): 

Vasym = (v 2 ym ~ 2v sym V k COS if/ + V^' 2 , (5) 

where v sym is given by Eq. (0, Vk = M co w/(M m + M 2 ), w is the 
kick velocity, and if/ is the angle between the kick velocity and 
the pre-explosion relative orbital velocity [note that Eq. (0 ne- 
glects the mass loss from the secondary star owing to stripping 
and ablation of the stellar material caused by the passing super- 
nova shock wave; see Tauris et al. (1 19991 ) for details]. 

It follows from Eq. (0 that v asym is maximum if the kick ve- 
locity attained by the stellar supernova remnant was oriented 
in the opposite direction to the relative orbital velocity, i.e. 
cos ip — - 1 . In this case 

Cym = <m + M co w/(M co + M 2 ). (6) 

Eq. shows that for HMXBs comprising a neutron star, i.e. for 
systems with M co /(M co + M 2 ) « 1, the effect of asymmetry 
of supernova explosion on the systemic velocity is modest and 
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could be substantial (larger than several tens of km s~') only for 
extremely large kick velocities, w > 1000 kmr 1 . It is gener- 
ally believed that neutron stars can obtain kick velocities of this 
magnitude. This belief got a support after Chatterjee et al. ( 120051 1 
derived a space velocity of ~ 1000 kms" 1 for the pulsar PSR 
B 1508+55 by measuring its proper motion and parallax. Still, 
the existence of a population of extremely high-velocity neutron 
stars can be explained without exploiting the idea of asymmet- 
ric supernova explosion - they simply could be the remnants of 
massive hypervelocity runaway stars (Gvaramadze 120071 120091 
Gvaramadze, Gualandris & Portegies Zwart 2008: 2009 1. 

Assuming that the neutron star in 4U 1907+09 attained a 
kick velocity of 200 - 400 km s" 1 at birth (typical of radio pul- 
sars; e.g. Hobbs et al. 120051 1 and adopting M co = 1.4 M Q and 
M 2 = 26 M G (Cox et al. 120051 ). one finds from Eq. © that the 
asymmetry of the supernova explosion can increase the systemic 
velocity only by ^ 10-20 kms" 1 (provided that the kick ve- 
locity was of favourable orientation; see above). We therefore 
conclude that the asymmetric supernova explosion did not con- 
tribute significantly to the space velocity of 4U 1907+09, unless 
the kick velocity was larger than 1000 km s" 1 . 

5.3. 4U 1907+09 and VelaX-1 

In many aspects, 4U 1907+09 is similar to VelaX-1: both sys- 
tems host OB supergiants of almost equal mass and radius, 
the orbital periods of both systems are almost equal to each 
other as well, and both systems are bow shock-producing run- 
aways. Below we compare the bow shocks associated with these 
HMXBs. 

The geometry of a bow shock generated by a supersonically 
moving star is characterized by two main parameters, the Mach 
number M = v*/c s (where v* is the velocity of the star relative 
to the ambient medium and c s is the sound speed of the ambi- 
ent medium) and the stand-off distance R s . The higher M (or the 
higher v„) the smaller the opening angle, a, of the bow shock, 
which is related to the Mach number through the relationship : 
a = 2 arcsin(l/M). In practice, however, this angle is difficult to 
determine owing to the low emission measure of the flanks of 
bow shocks. Still, a comparison of Fig.Q~]and Fig |2] suggests that 
the opening angle of the bow shock generated by 4U 1907+09 is 
much smaller than that of the bow shock associated with Vela X- 
1, meaning that the peculiar velocity of the former system is 
much larger than that of the latter one (see also below). 

Let us estimate the space velocity of VelaX-1 using the 
proper motion of this HMXB from the new reduction of the 
Hipparcos data by van Leeuwen ( 120071 1. p a co&8 = -5.48 ± 
0.35 masyr" 1 ,/^ = 8.79 + 0.40 masyr" 1 , and the generally ac- 
cepted distance to VelaX-1 of 1.9 kpc (Sadakane et al. 119851 ). 
Using the Galactic constants and the solar peculiar motion 
adopted in Sect. 15. 11 we found the transverse peculiar velocity 
of VelaX-1: 
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-35.4 ± 3.4 kms" 



23.8 ±3.3 kins" 



(7) 



or v tt = 42 7 + 3.4 kms" 1 (cf. Moffat et al. [19981 Huthoff & 
Kaper I20021 ). To this velocity one should add the peculiar ra- 
dial velocity v ra d = -27 kms" 1 , inferred from the heliocentric 
one of -4 km s" 1 (Gies & Bolton [19861 ), so that the total (three- 
dimensional) velocity of VelaX-1 v» - 50 kms -1 , and its vector 
is inclined to the plane of the sky by an angle of 32°. 

The second parameter characterizing a bow shock is the 
stand-off distance, i.e. the minimum distance from the star at 
which the ram pressure of the stellar wind is balanced by the 



ram pressure of the ambient interstellar medium (e.g., Baranov, 
Krasnobaev & Kulikovskii l 19711) : 



Rs = 



Anpvl 



1/2 



(8) 



where M and v w are the stellar mass-loss rate and the wind ve- 
locity, p = 1 Am^n, inn is the mass of a hydrogen atom, and n is 
the number density of the ambient medium. R s is related to the 
observed angular separation between the apex of the bow shock 
and the star, 9, through the relationship: 



G=/3y(i)Rs/d, 



(9) 



where (5 > 1 is a factor of order unity taking into account the 
finite thickness of the layers occupied by the shocked gas coming 
from the ambient medium and the stellar wind [ see, e.g., Fig. 1 
in Comeron & Kaper (119981 for schematic of the multi-layer 
structure of a wind bow shock], 



r(0= j_j 3 (i--4 

sin; V V tan;/ 



(10) 



is a geometric factor taking into account the inclination of the 
bow shock to the plane of the sky (see Wilkin 119961 ). where 
i = arctan(v ra dMr) i s me angle between the vector of the space 
velocity of the star and the plane of the sky. Note that for small 
inclination angles (; < 45°) the factor y is ^ 1. 

Using Eqs (fHTt-dTOb. one can estimate the number density of 
the ambient medium: 



n 0.5 cm" 3 (J3yY 



M 



lO^Moyr- 1 
x 



( *2 ) 

1 2000 kms- 1 / 



^obs 

Tpcj 1 100 kms" 1 



(11) 



where R a ^ = 9d. Inserting parameters of 4U 1907+09 and 
VelaX-1 (given in Table[TJ in Eq. (fTTb . one finds the number 
density of the interstellar medium surrounding these HMXBs, 
=s 2 and lOcnT 3 , respectively (we assumed here that /3y ~ 1). 
Note that despite the much larger separation of VelaX-1 from 
the Galactic plane 130 pc versus 30 pc for 4U 1907+09) it 
is located in a denser medium than 4U 1907+09, implying that 
VelaX-1 met a density enhancement (a cloudlet) on its way (cf. 
Gvaramadze & Bomans 2008 1 Gvaramadze et al. 2009 1. This im- 
plication is supported by the presence of an extended (of angular 
radius of ~ 0°2) IRAS source at the position of VelaX-1, sug- 
gesting that this source might be a Stromgren zone created by the 
ultraviolet emission of VelaX-1 (Kaper et al. 119971 ). Assuming 
that the gas in the H n region is fully ionized and its temperature 
is 10 4 K, one can obtain an estimate of the gas number density 
(e.g. Lequeux l2005l) 



n =: 15 cm 



7pc 



-3/2 



S(0) 



3 x 10 48 photons s~ 



1/2 



(12) 



where R$ t is the Stromgren radius (derived from the angular ra- 
dius of the Hn region and the distance to VelaX-1 of 1.9 kpc) 
and S(Q) is the total ionizing-photon luminosity of an B0. 5 lab 
star (cf. Martins, Schaerer & Hillier 12005) . The good agreement 
between the two independent estimates derived from Eqs dTTT > 
and dT2b supports the interpretation of the extended infrared 
emission around VelaX-1 and its bow shock as the Hn region 
[Kaper et al. 119971 see also Fig. 5 in Gvaramadze & Bomans 
( 2008 ) for another example of a bow shock within an H n re- 
gion]. 
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Table 1. Summary of parameters related to 4U 1907+09 and VelaX-1. 



HMXB 


Spectral type 


d (kpc) 


v»(kms-') M(10^M o yr') 


v w ( km s 1 ) 


fiobs (pc) 


n (cm 3 ) 


4U 1907+09 


09.5 lab" 


4 


160 6" 


1690" 


0.8 


2 


VelaX-1 


B0.5 Iab c 


1.9 


50 2" 


1100" 


0.5 


10 



Notes. (fl) Nespoli et al. l2008l <*> Cox et al. l2005l (c) Conti [19781 (rf) Watanabe et al. l2006l (c) Prinja, Barlow & Howarth [T990l 



It should be noted that the geometric factor y(i) [given by 
Eq. (TTOb l was derived in the thin-shell approximation, i.e. under 
the assumption that the thickness of the region occupied by the 
shocked gas is << R s (Wilkin [19961 >. This approximation, how- 
ever, could be valid only for bow-shocks generated in the very 
dense medium (e.g. cometary ultra-compact H n regions; Mac 
Low et al. l 1 99 11 > or for bow shocks produced by stars with dense, 
slow winds (e.g. red supergiants). In both cases, the shocked gas 
could be dense enough to rapidly cool and collapse into a thin 
shell. The thin-shell approximation hardly can be used to de- 
scribe the shape of adiabatic bow shocks, i.e. the bow shocks 
generated by massive stars possessing fast, hot winds and/or 
moving through the space with high peculiar velocities. For 
these bow shocks, the thickness of the layers occupied by the 
shocked stellar wind (the inner layer) and the shocked ambient 
medium (the outer layer) constitutes a significant fraction of R s 
(Raga et al. [19971 Comeron & Kaper [T99^b . 

Because 4U 1907+09 and VelaX-1 both contain supergiant 
stars with strong winds (see Table[TJ, one can expect that at least 
the inner layer of the bow shocks generated by these HMXBs 
is adiabatic. To check this, we use the following back-of-the- 
envelope estimates. The shocked gas in a bow shock will col- 
lapse into a thin shell if the characteristic time-scale for radiative 
cooling of the gas, 



^cool ~ 



3kT 
n ps A(T) ' 



(13) 



is much shorter than the characteristic dynamical time-scale of 
the bow shock, 



t ^ 

'dyn 



V* 



(14) 



where k is the Boltzmann constant, r ps and « ps are the initial 
post-shock temperature and number density, and 



A(T) = 



(0,T < 10 4 K 
1.0 x 10- 24 r (I55 ,10 4 K< T < 10 5 K 
6.2 x 10- I9 r- (X6 , 10 5 K < T < 4 x 10 7 K 
2.5 x 10" 27 r 05 ,r >4x 10 7 K 



(15) 



is the cooling function (Raymond, Cox & Smith 119761 
Cowie, McKee & Ostriker I198U . For the (inner) layer of 
the bow shock occupied by the shocked stellar wind, T ps 
1.4 x 10 7 (v w /1000 kms~') 2 and « ps = 4« w , where « w 
M/(5.67toth/? 2 v w ). Using Eqs. (TT3T> - (fT~5T> and parameters of the 
stellar winds from TableQ] one finds that f coo i >> fdyn for both 
bow shocks, i.e. the inner layer of the bow shocks is adia- 
batic. Similarly, adopting for the shocked ambient gas T ps 
1.4xl0 5 (v„/100 kms" 1 ) 2 and« ps = An, one finds that f coo ] ~ fd yn 
for 4U 1907+09 and t coo \ « f dyn for VelaX-1, so that the outer 
layer of the bow shock produced by 4U 1907+09 is adiabatic, 
while that of the bow shock around VelaX-1 is radiative. From 
this it follows that the shape of both bow shocks cannot be 
described by the thin-shell solution of Wilkin (1996). Indeed, 



the Wilkin's solution predicts that the radial extent of a bow 
shock in the direction perpendicular to the direction of motion 
R(tt/2) = y(n/2)R s =* 1.7J? S , while from Fig.[TJand Fig.|2]we 
found R(n/2) l.5R s for 4U 1907+09 and R{n/2) = 2AR S for 
VelaX-1. 

It is also instructive to compare the geometry of the bow 
shocks produced by 4U 1907+09 and Ve laX-1 with that of 
model bow shocks from Comeron & Kaper ( 1998), namely with 
their case D and case E bow shocks (generated by stars with 
space velocities of 50 and 150 kins -1 , respectively, which are 
almost identical to those of Vela X-l and 4U 1907+09). One can 
show that like in the case of 4U 1907+09 and Vela X-l, the in- 
ner layer of both model bow shocks is adiabatic (cf. Comeron 
& Kaper 1998), while the outer one is radiative in the case D 
bow shock and adiabatic in the case E bow shock. It is not sur- 
prising, therefore, that the bow shock around Vela X-l shows a 
close similarity to the case D bow shock, while the shape of the 
bow shock associated with 4U 1907+09 reminds that of the case 
E bow shock [cf. Fig.|2]and Fig.[T]with the upper and the lower 
panels of Fig. 13 in Comeron & Kaper (1998), respectively]. 
The more thorough study of the bow shocks associated with 
4U 1907+09 and VelaX-1 requires numerical simulations [sim- 
ilar to those carried out by Comeron & Kaper (JT998 »], which is 
beyond the scope of the present paper. 

The inference that the outer layer of the bow shock around 
Vela X-l is radiative is consistent with the observational fact 
that just this of the two bow shocks manifests itself in Ha. On 
the other hand, the non-detection of the Ha emission from the 
bow shock associated with 4U 1907+09 could be used to some- 
what constrain the distance to this HMXB. The space veloc- 
ity of 4U 1907+09 approximately scales with the distance as 
~ d, so that n ~ d 4 , and correspondingly, t coo \ oc d (for 
10 5 K < T < 4 x 10 7 K) and t dyn oc d°. The strong dependence 
of fcooi ° n d implies that the outer layer of the bow shock would 
be radiative, while the bow shock itself would be visible in Ha 
if the distance to 4U 1907+09 is shorter than 4 kpc, say compa- 
rable to that of Vela X-l (^ 2 kpc). This implication, however, 
suffers from the large error bar in the proper motion measure- 
ments, which causes the correspondingly large uncertainties in 
estimates of the space velocity of the HMXB and the cooling 
time-scale of the bow shock. And on the related note, for dis- 
tances much larger than 4 kpc, the outer layer of the bow shock 
would remain adiabatic, while the space velocity of the sys- 
tem would become too high to be explained within the binary- 
supernova scenario (see Sect. 15. 2\ . 

Finally, we note that the IRAC 8 fim image of the bow shock 
produced by Vela X-l (right panel of Fig.|2j shows filamen- 
tary (cirrus-like) structures beyond the main body of the bow 
shock. We speculate that these filaments are due to interstellar 
dust grains aligned with the local interstellar magnetic field and 
heated by the radiation of Vela X- 1 . 
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6. Summary and conclusion 

We have searched for bow shocks around 5 HMXBs (GX 301-2, 
4U 1 907+09, CygX-1, Vela X-l and V615Cas) using the Spitzer 
Space Telescope archival data. We found two bow shocks, one 
of which (produced by Vela X-l) was already known, while the 
second one (generated by 4U 1907+09) was detected for the 
first time. The detection of the bow shock around 4U 1907+09 
provides strong evidence that this HMXB is a runaway sys- 
tem, supporting the general belief that most of HMXBs pos- 
sess high space velocities caused by supernova explosion of 
one of the components of the progenitor massive binary. To 
prove the runaway nature of 4U 1907+09, we have measured its 
proper motion using archival data over a total time baseline of 
50 yr. Although suffering from large error bars, our measure- 
ments support the runaway interpretation of 4U 1907+09, and 
suggest that this system was expelled from the Galactic plane 
about 10 5 yr ago. Using the observed parameters of 4U 1907+09, 
we found that the high space velocity of this HMXB could be 
well explained within the framework of the symmetric binary- 
supernova explosion. We found also that the possible asymme- 
try in the supernova explosion does not contribute significantly 
to the space velocity of 4U 1907+09. We have compared the 
bow shocks produced by 4U 1907+09 and Vela X-l and showed 
that the cooling of the shocked stellar wind is inefficient in both 
shocks. We showed also that the layer of the shocked ambi- 
ent gas is adiabatic in the former bow shock and is radiative in 
the latter one. This allowed us to explain qualitatively why the 
shapes of both bow shocks do not follow the thin-shell solution 
by Wilkin < fl996b . 

To conclude, we expect that the future proper motion and 
parallax measurements for 4U 1907+09 with the space astrom- 
etry mission Gaia will considerably improve the estimate of its 
space velocity, and thereby will allow us to better understand 
the evolution and pre-supernova configuration of its progenitor 
binary system. With the new proper motion and distance mea- 
surements in hands we will also be able to refine the parameters 
of the bow shock around 4U 1907+09 and to confront them with 
numerical simulations. 
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